Negative oxygen vacancies in Hf02 as charge traps in high-k stacks 
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We calculated the optical excitation and thermal ionization energies of oxygen vacancies in m- 
Hf02 using atomic basis sets, a non-local density functional and periodic supercell. The thermal 
ionization energies of negatively charged V~ and centres are consistent with values obtained 
by the electrical measurements. The results suggest that negative oxygen vacancies are the likely 
candidates for intrinsic electron traps in the hafnum-based gate stack devices. 



Hafnium based oxides are currently considered as a 
practical solution satisfying stringent criteria for integra- 
tion of high-k materials in the devices in future technol- 
ogy nodes. However, high-k transistor performance is 
often affected by high and unstable threshold potential-, 
Vt, and low carrier mobilitj*^. These effects arc usually 
attributed to a high concentration of charge traps and 
scattering centers in the bulk of the dielectric and/or at 
its interface with the silicon channel. Although the re- 
ported trap densities vary greatly with fabrication tech- 
niques, the majority of data point to existence of a spe- 
cific intrinsic shallow electron centre common to all Hf02 
based stacks while some extrinsic defects, such as Zr sub- 
stitution in Hf02, have also been considered^. 

Oxygen vacancies are dominating intrinsic defects in 
the hulk of many transition metal oxides including Hf02 
and Zr02, and are thought to be also present in high con- 
centrations in thin films. However, in spite of numerous 
experimental studies, evidence relating oxygen vacancies 
to measured characteristics of interface traps in high-k 
stacks is still mostly circumstantial. Therefore, accurate 
theoretical characterization of these defects is highly de- 
sirable. 

Previous theoretical calculations of oxygen vacancy in 
Hf02 and Zr02 reported the ground state properties ob- 
tained within local or semilocal approximations to den- 
sity functional theory (DFT) methods (see^ for a re- 
view). This approach, however, significantly underesti- 
mates band gaps, which hampers determining energies of 
defect levels with respect to the band edges and precludes 
identifying shallow defect states^iSii. As a result, most of 
the early local DFT calculations (except, perhaps, refiS) 
failed to predict unambiguosly negative charge states of 
oxygen vacancy in Hf02. Significant improvement was 
achieved by Robertson et almi^^ who used screened 
exchange approximation to predict vacancy energy lev- 
els including V~ charge state. However, these calcu- 
lations were performed using a small periodic supercell 
and therefore corresponded to extremely high vacancy 
concentrations. The quality of the functional used is also 
largely unknown and needs independent verification. 

In this work we used much bigger supercells and a non- 
local functional to calculate optical excitation and ther- 



mal ionization energies of oxygen vacancies in five charge 
states. To relate these energies to experimental data we 
distinguish optical absorption/refiection type measure- 
ments involving Frank- Condon type (vertical) excita- 
tions, and electrical thermal de-trapping measurements, 
where phonon-assisted electron excitations are accompa- 
nied by strong lattice relaxation. We focus on the results 
of so called de-trapping electical measurementsiiiiSii^ 
which are interpreted in terms of thermal ionization of 
shallow electron traps and demonstrate that the inter- 
pretation is consistent with thermal ionization of nega- 
tive V~ and vacancies. 

Our periodic non-local density functional calculations 
were carried out using atomic basis set and a B3LYP 
hybrid density functionali4 implemented in the CRYS- 
TAL03 codoi^. This method reproduces band gaps of a 
range of transition metal oxidesi^ and allows us to carry 
out geometry optimization of defect structures, calcu- 
late defect electronic properties and excitation energies 
within the same method. We used all electron basis- 
set for oxygen atomaii, and an s,p,d valence basis set 
for Hf in conjunction with the relativistic small core ef- 
fective potential due to Stevens et alm^. All calculations 
were carried out in the 96 atoms supercell of a monoclinic 
(m)-Hf02 with a mesh of 36 k points in the irreducible 
Brillouin zone. The defect structures were optimized to 
the atomic forces below 0.03 cV A~^. The compensat- 
ing uniform background potential method was used for 
charged defects calculationaSiiS. The positions of occu- 
pied and un-occupied single-electron defect levels in the 
gap of m-Hf02 predicted by periodic B3LYP calculations 
are presented in Figure ^ The results of the calculations 
can be summarized as follows; 

1. The calculated single-electron band gap in m-Hf02 is 
6.1 eV. We note, that single-electron estimate neglects 
excitonic effects, and effects of phonon broadening which 
reduce measured optical band gap of Hf02 by 0.8 eV—. 
Taking this into account, predicted Eg is in good agree- 
ment with the experimental optical gap of 5.6-5.9 eV—. 

2. Oxygen vacancy in the monoclinic Hf02 may exist in 
five charge states +2, +1, 0, -1, -2 corresponding to up 
to four extra electrons in the vicinity of the vacant O^^ 
site. 
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FIG. 1: Defect energy diagram for the 4-coordinated oxygen 
vacancy in the m-Hf02. 



3. The preferential site for the oxygen vacancy formation 
in m-Hf02 depends on the charge state: V^^ and V+ are 
more stable in the 3-fold coordinated sites, whereas the 
V°, V~, and are energetically more favorable at the 
4- fold coordinated sites (see also refsi^^S^). The difference 
in formation energies between 3- and 4-coordinatcd sites 
of negatively charged vacancies is about 0.2 cV. 

4. The density of electrons in the +1 and zero charge 
vacancy is peaked at the vacant site with significant ad- 
mixture of d-states of the nearest neighbor hafnium ions, 
similar to plane wave DFT calculation a^i^^ . The single- 
electron states of these electrons arc located roughly in 
the middle of the gap (Fig. 1), significantly lower than 
those reported by Robcrtson^i2iifi 

5. Both one and two extra electrons added to a neutral 
vacancy form more diffuse asymmetric states localized 
mainly on 3 (out of 4) Hf ions surrounding the vacancy. 
The ground state of the is spin singlet, which is 0.2 
eV lower then the spin triplet configuration (cf. Fig. 

6. The change of the vacancy charge state is accompa- 
nied by significant displacements of the nearest neighbor 
(NN) Hf atoms and next shell of oxygen atoms (NNN) 
(see also refi^. The NN Hf ions displace from their per- 
fect lattice position approximately symmetrically. This 
displacement is away from the and V+ by 11% and 
5% of a typical Hf-0 distance, but towards the V~ and 
V^^ by 4% and 8% ,respectively. The displacements of 
the NNN oxygen ions are substantially smaller and di- 
rected away from the negative vacancies. 

7. The character of electron density distribution, strong 
lattice relaxation and relatively shallow single-electron 
levels (Fig. suggest that trapping of extra electrons in 
negative V~ and vacancies is essentially polaronic 
in nature. The third and forth electrons induce strong 
lattice polarisation, which in turn, creates the potential 
well for these electrons. Now we discuss a relation be- 
tween the single-electron spectrum of Fig. ^ and exper- 
imental trap ionization measurements. For optical type 
measurements, such as optical absorption, and reflection 
ellipsometrjiS^, using energy differences between single- 
electron levels can provide reasonable accuracy^iS^. The 
situation is different in electrical measurements, which 
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FIG. 2: Schematic potential energy surface (PES) for optical 
and thermal detrapping processes. Generalized displacement 
of the Hf atoms, Q, and energy scale is chosen to represent 
ionization of V^^ vacancy. The band diagram schematic illus- 
trates optical and thermal ionization processes, respectively 
without (B*) and with (B) lattice relaxation. 

include thermal excitation of trapped electrons into the 
conduction band of HfO^ii (see inset of Fig. EJ- 

The difference between the optical and thermal ion- 
ization can be qualitatively explained by the potential 
energy surface diagram depicted in Fig. |21 and earlier 
discussed for other systems (e.g. in refii). Fast optical 
ionization of a trap of charge q (state A) into the bottom 
of the conduction band (state B*) corresponds to a Frank- 
Condon type transition, Eopt, when no lattice relaxation 
associated with the removed charge occurs. In contrast, 
a much slower thermal ionization of a trap, E^;,, is a 
phonon assisted process, during which the system trans- 
fers into a fully relaxed trap state q + 1 and the electron 
delocalized at the bottom of the conduction band (state 
B). The ionized trap induces lattice polarization, which 
in case of the oxygen vacancy is mainly associated with 
strong displacements of the NN Hf ions. This generalized 
relaxation coordinate is denoted as Q in the diagram in 
Fig. 121 Et/i is approximately given by: 

Eth — Eopt — Erel, (1) 

where E^e; is the lattice relaxation energy. Optical and 
thermal ionization energies can be calculated using the 
total energies of systems in different charge states as fol- 
lows: 

EoptiV) = Eg{V'^+^) - E.iV) + E-- E"; (2) 

EreliV'^+') = E,+ ,iV^+') - E,{V^+^). (3) 

The notations in Eqs. I2I3I have the following meaning: 
E° - the total energy of the perfect Hf02 crystal; E~ - 
the total energy of the perfect Hf02 crystal with an elec- 
tron at the bottom of the conduction band; Eq(V'') - the 
total energy of Hf02 with the vacancy in a charge state q 
(q= -1-2, -1-1,0,-1,-2) in the optimized geometry; Eq(V''+^) 
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TABLE I: Optical excitation energies, Eopt, relaxation ener- 
gies. Ere!, and thermal activation Energies, Eth, for oxygen 
vacancies in m-Hf02 calculated according to Eqs. I1I3I 



q 


Eopt Erel Eth 


v+ 
v« 

V- 


3.33 1.01 2.32 
3.13 0.80 2.33 
1.24 0.48 0.76 
0.99 0.43 0.56 



- the total energy of the vacancy in the charge state q+1 
but at the equilibrium geometry corresponding to the va- 
cancy in the charge state q. We assume that the vacancy 
is well separated from the interface, and the optical and 
relaxation energies discussed here, depend only on the 
bulk properties of Hf02 and not on the Si or metal band 
alignment. The calculated values of Eopt, ^rei and l^th 
are summarized in Tabled We note that optical ioniza- 
tion energies, Eopt calculated as total energy differences 
in Eq. |21 are close to single-electron energy differences 
(Fig. This is to be expected, since the total energy 
differences for the systems with N and N±l electrons in 
the DFT calculations are related to the single electron 
energies of the highest occupied and lowest unoccupied 
statesSi. 

The thermal ionization energies Eth are, however, 0.5 



to 1.0 eV smaller than the optical energies due to the 
large lattice relaxation associated with the change of the 
charge state of the vacancy. Although the single parti- 
cle energies of V~ and are very similar (Fig. 
thermal ionization energies of these defects differ by 0.2 
eV. These values arc consistent with 0.35, 0.5 eV activa- 
tion energies extracted from thermal de-trapping kinet- 
ics measurementsiiiSii^. On the other hand, very large 
de-trapping energies for the V"*" and V° vacancies rule 
out these species as possible shallow electron traps. We 
should note that accounting for the thermal broadening 
of the defect levels and band tails, which will be discussed 
elsewhere, would improve the agreement between the cal- 
culated thermal ionization energies and the experimental 
energies even further. 

To summarize, we have calculated the positions of 
single-electron levels for five charge states of oxygen 
vacancy in m-Hf02 and related them to the existing 
experimental data. The results of trap discharging 
measurementsiiii^ii^ are consistent with thermal ioniza- 
tion of negatively charged V~ and oxygen vacancies. 
These results further support the common assumption 
that oxygen vacancies are likely candidates for intrinsic 
electron traps in these devices and suggest that negative 
oxygen vacancies can be responsible for Vt instability. 
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